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Of the available classes of insulators which have been shown to contain topologically non-trivial 
properties one of the most important is class All, which contains systems that possess time-reversal 
symmetry T with T 2 — —1. This class has been the subject of significant attention as it encompasses 
non-trivial Z2 topological insulators such as the quantum spin Hall (QSH) state and the 3D strong 
topological insulator. One of the defining properties of this system is the robustness of the state 
under the addition of disorder that preserves T. In this letter, we explore the phase diagram of 
the disordered QSH state as a function of disorder strength and chemical potential by examining 
the entanglement spectrum for disordered class All symplectic systems. As for the case of the T 
breaking Chern insulator we show that there is a correspondence between the level-spacing statistics 
of the Hamiltonian and that of the level spacing statistics of the entanglement spectrum. We observe 
a feature in the statistics of the entanglement spectrum that aids the identification of delocalized 
states and consequently critical energies across which phase transitions occur. 



The classification of topological insulators and super- 
conductors via the periodic table |TH3] has garnered much 
interest over the last few years. One of the most ex- 
citing results to have come out of the classification is 
the connection between the ten symmetry classes of 
gapped free-fermion systems and the ten random ma- 
trix ensembles studied by Altland and Zirnbauer[2, 4j. 
The idea that each symmetry class of topological states 
has a corresponding random-matrix ensemble enables 
one to carefully study the properties of disordered 
topological insulators, and their corresponding localiza- 
tion/delocalization transitions between topological and 
trivial insulators. As an example, the phase diagram of 
the disordered Chern insulator [5 , which belongs to the 
Wigner-Dyson unitary ensemble denoted by class A, has 
been calculated in Ref. [SJ The other experimentally rele- 
vant class is All which is the symplectic class which con- 
tains insulators with time-reversal symmetry T such that 
T 2 = -1. This class is known to support both 2D[?HS] 
and 3D[10 topological insulators classified by a Z2 topo- 
logical invariant. 

Our focus will be on the 2D insulator of class All which 
is also known as the quantum spin Hall (QSH) state. The 
properties of the disordered QSH state have been previ- 
ously studied in Ref. by numerically evaluating the 
Z2 topological invariant while varying disorder strength, 
chemical potential, and the nominal insulating gap. We 
seek instead to study this system via the entanglement 
spectrum O H2TTH] in the hope of developing new tools 
to study phase transitions in disordered systems by only 
considering ground state properties and not the excita- 
tion spectrum. Recently it has been shown that many 
properties of the entanglement spectrum mimicked prop- 
erties of the real energy spectrum; specifically the statis- 



tics of nearest-neighbor level spacings[5]. There are two 
additional parameters in the calculation of the entangle- 
ment spectrum with respect to the Hamiltonian spec- 
trum: (i) the choice of the ground-state (single-particle 
energy level filling) (ii) the choice of the Hilbert space 
partition/cut. In this Letter we will show that, in the 
context of topological insulators, studying how the statis- 
tics of the entanglement spectrum varies when the single- 
particle filling varies can reveal regions of delocalized 
states in the energy spectrum in the presence of quenched 
disorder. 

We will consider an effective square-lattice model for 
the 2D QSH Hamiltonian[9 as a representative Hamil- 
tonian for the All symmetry class. The Fourier trans- 
formed clean Hamiltonian is simply 

= £4 (T *•<-*> m 

with 

H(k) = d(k) ■ t (2) 

where = (cyt-j- cb^ ca\. cbi) T with A and B repre- 
senting the orbital degree of freedom, f = (t x ,t v ,t z ) 
are the Pauli matrices for the orbitals, and d(k) — 
(sinfc^, smk y , 2 — m — cosk x — cosk y ). The topological 
state of interest occurs when < m < 2, and, for exam- 
ple, a trivial state occurs when m < 0. A model similar to 
this was used as an effective model for HgTe / CdTe quan- 
tum wells [9 . While this model represents a time-reversal 
invariant topological insulator, its random matrix prop- 
erties when time-reversal invariant disorder is added ex- 
actly match that of two copies of the Chern insulator 
and thus falls into class A. To study the properties of 
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the symplectic class All we need to deform this model, 
which we will do by adding spin-mixing terms which pre- 
serve T, notably a term which is due to the bulk inversion 
asymmetry (BIA) of the HgTe/CdTe zinc-blende lattice 
structure. We will choose the BIA term 
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where A is the strength of the BIA which was derived 
in Ref. 1151 This term mixes spin-up and spin-down, 
preserves T, but breaks inversion symmetry; in turn it 
deforms our model away from a pathological case so that 
it represents class AIL The last piece of the Hamiltonian 
is the disorder which we take to be on-site potential dis- 
order with a random strength uniformly distributed over 
[-W/2, W/2] where W is a real parameter. 

With the Hamiltonian established, we now examine 
the Hamiltonian energy level statistics and the resultant 
entanglement spectra statistics. In collecting statistical 
data, we use a lattice model with size 40 x 40 with peri- 
odic boundary conditions imposed on both the x and y 
directions so that there are no low-energy edge states and 
the system is completely gapped (in the clean limit) . We 
use exact diagonalization to calculate the energy spectra 
and the entanglement spectra via Peschel's methodj? ]; 
essentially calculating the two point correlation function 



Ci 



(etc, 



(4) 



as a function of the lattice positions i, j where the brack- 
ets ( ) mean the ground-state expectation value for a 
single disorder configuration. If one restricts i,j to lie 
in a subset A of the entire lattice then one can extract 
the single-particle entanglement spectrum by diagonaliz- 
ing this truncated matrix. For our calculations we took 
the region A to be a rectangular subset of size 20 x 20. 

Let us first examine the nearest-neighbor level-spacing 
statistics for the energy spectrum|16|. To calculate the 
distribution of level spacings one picks an energy Eq in 
the spectrum and calculates the nearest-neighbor level 
spacings of states in a small window around Eq and then 
combines the level spacing data for many different ran- 
dom disorder configurations to form a distribution for 
each choice of Eq. In this work, we use a window which 
considers states 3 levels above and below Eq. In Fig. 
[TJa) we plot the normalized variance ((s 2 ) — (s) 2 ) of this 
distribution for a large range of choices for Eq covering 
the entire spectrum of the Hamiltonian. This data was 
taken for W — 5 and A = 0.3 and averaged over 1000 
disorder configurations. Our model is particle-hole sym- 
metric (but only on average disorder is included) and 
the resultant distribution variances are also symmetric. 
For class All we expect two different types of signa- 
tures: (i) a Poisson distribution (Pp = e~ s ) for states 
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Figure 1. (a) Plot of the level statistics of a QSH insulator 
with a BIA term of A = 0.3 and a random disorder potential 
of magnitude W = 5. In plots (b)-(d) we show the resultant 
entanglement spectrum calculated at chemical potentials of 
(b) fi = 3 (c) fi = 2 (d) fj, = 1 
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Figure 2. Plot of the level statistics of a QSH insulator with 
W = 1 for various values of the BIA term. (a)A = 0.0 
(b)A = 0.0 (c)A = 0.0 (d)A = 0.0 We clearly see that as A in- 
creases the width of the delocalized region of states increases 
proportionally. 



which are localized with a variance which approaches 1 
in finite-size (ii) a symplectic Wigner-Surmise distribu- 
tion (Pgse = Ns A cxp f — i with N a normalization) 
with a variance of which approaches ~ 0.1045. We see ev- 
idence for both types of distributions depending on the 
value of the chemical potential. Near the band edges, we 
see the variance trending toward 1; in fact, if the system 



size is increased the variance away from the center of the 
spectrum more quickly asymptotes to a value of 1 indi- 
cating mostly localized states. However, there are two 
regions (near E — ±1) which have variances very close 
to the Wigner-Surmise value (red- vertical line). Com- 
pared to similar data of two of the authors shown in Ref. 
|5] we notice a difference in that there seems to be a rel- 
atively wide band of states that reach variance ~ 0.1045 
compared to a much narrower range for the Chern insu- 
lator. We take this to be a signature of a critical metallic 
region (as opposed to a critical point) separating the triv- 
ial from topological insulator phases. For the disordered 
All class such a region is expected [TT]. As an aside let us 
note that if we tune A — > the delocalized region thins, 
and if we increase A (as shown in Fig. [2] for W = 1) the 
width increases which shows that the band of delocal- 
ized states is due to the symplectic nature of the random 
matrices. Thus for class All one must pass through an 
intermediate metallic/critical phase when going from the 
disordered topological insulator phase to the disordered 
trivial insulator phase via the levitation and annihilation 
transition mechanism. 

Now this behavior is to be compared with the resultant 
entanglement spectra of the same system with ground- 
state fillings given by several chemical potentials shown 
in Figs [ljb)-(d). Here we see that the EnS also has a 
regions with large variances of order 1 and regions with 
variances that closely approach 0.1045. In Fig. [ljb) the 
filling at fj, = 3 is such that the states at the Fermi- 
level are localized, and the system is in a trivial insulator 
phase. We see that the variance of the entanglement en- 
ergy levels (denoted £j) derived from the eigenvalues of 
the truncated Cu (denoted Q) result in a large variance 
and exhibit a Poisson-like distribution for essentially ev- 
ery choice of the entanglement energy £q. At fi = 2 the 
picture is less clear but it appears that some regions of 
entanglement energy are approaching the delocalized dis- 
tribution, though the curves do not actually touch the 
asymptotic value. This does not just seem to be a finite- 
size effect since as we approach closer to the delocalized 
band of states the variances actually get much closer to 
the asymptotic Wigner-Surmise value. For /i = 1 where 
the filling is very nearly in the band of delocalized states 
we see the variance of the entanglement energy spacings 
flattening out to the delocalized value. It should be noted 
that a similar effect was seen in the case of the Chern 
insulator where, when tuned to into the band of delo- 
calized states the variance of the entanglement energies 
completely flattened out to the Wigner-surmise value [TT], 
It appears that this phenomenon is a general property 
(at least of topological insulators) as we tune the chemi- 
cal potential through delocalized states. We have verified 
that this behavior with larger system sizes of 50 x 50 with 
a corresponding 25 x 25 region A in order to rule out the 
possibility of significant finite size effects skewing the re- 
sults of the entanglement. We also note that all of the 
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Figure 3. Plots of the variance of the level spacings for the full 
spectrum corresponding to random disorder strengths of (a) 
W=5.1, (b) W=5.625, (c) W=6 and (d) W=W. We see that 
as disorder strength is increased we pass the quantum critical 
point and the system passes from topological to metallic then 
to insulating. In (e)-(h), we plot the entanglement spectra 
for /i=0 corresponding to the disorder strengths considered in 
(a)-(d) with a BIA term of A=0.3. 



entanglement variances have a similar shape with a peak 
near £ = 0. This property also occurred for the Chern in- 
sulator and is not completely understood. Our intuition 
indicates that this peak is arising from "local entangle- 
ment." States that give rise to the low energy part of 
the entanglement spectrum (i.e. those states which are 
highly entangled) arise from two mechanisms: (i) they 
are delocalized states which have significant probability 
to be in region A and its compliment (ii) they are local- 
ized states which, by chance, are localized near the cut 
and contribute to the entanglement in spite of their local 
nature. A method to remove local entanglement from 
these calculations is still lacking, but in principle should 
make these calculations easier to interpret. 

Having considered the properties of the entanglement 
statistics as a function of chemical potential we now wish 
to examine the dual case of a fixed chemical potential but 
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varying disorder strength (still at A = 0.3). To calculate 
the entanglement spectrum we will fix our ground-state 
filling at /i = and focus near the region of disorder 
strength where the system exhibits phase transitions be- 
tween a topological insulator, metal, and trivial insulator 
at this energy. For our model this occurs when the dis- 
order strength W satisfies 5 < W < 6. We show energy 
and corresponding entanglement variances (for a filling 
of n = 0) in Fig. [| For Figs, ^a), (b), (c) and (d), we 
see the variances for the energy spectrum in the topolog- 
ical insulator (W = 5.1), metal (W = 5.625), and trivial 
insulator (W = 6 and W— 10) phases respectively. We es- 
timate the critical point separating the metal phase from 
the trivial insulator occurs around W = 5.85 which is the 
point where the energy variance curve near E — lifts 
off the delocalized value. The phase transition occurs as 
the two bands of delocalized states merge near E = 
and eventually annihilate in a manner consistent with 
the scenario of levitation and pair annihilation [17]. As 
these are clearly bands of delocalized states there is a fi- 
nite range of W over which these states are annihilating 
which gives rise to the critical metallic region. We see 
clear features in the corresponding entanglement spec- 
tra, with the most interesting one seen in Fig. [3jf ) which 
shows a flat variance curve for the entanglement energies 
in the case where the chemical potential intersects the 
critical metallic region at E = 0. The variances for the 
two insulating phases in Figs. [3je) (topological) and (g) 
(trivial) may be distinguished by noting that the topo- 
logical phase contains regions of extended states touching 
the prescribed variance while the trivial phase exhibits 
bands of states which have a variance which sits above the 
requisite variance. This is exacerbated in Fig. (3jh) where 
the entanglement spectrum sits well above the variance 
denoting the presence of delocalized states. Figs. [3jd) 
and (f) show evidence of contributions from local en- 
tanglement which produces the central maximum. This 
maximum seems to be a feature of the gapped phases as 
it does not appear when the system is tuned to a critical 
metallic region either by varying chemical potential, or in 
this case, varying disorder strength until a critical region 
hits the fixed /x. 

As hinted in Ref. [5], and shown to be a more gen- 
eral feature in this work, the flattening of the variance of 
the entanglement energy spacing is an indicator that the 
filling level is tuned to a region of delocalized states and 
can be used to identify critical regions effectively in order 
to map out a phase diagram. Using this technique, we 
can glean information on the resultant topological phase 
transitions from a careful analysis of the EnS. We are now 
in a position to discuss the phase diagram of our disor- 
dered QSH model. In Fig. [4] we plot the phase diagram 
as a function of random disorder strength and chemical 
potential for A = 0.3. We see a connected region of the 
quantum spin Hall phase separated from an outer region 
of trivial insulator by an intervening metallic region. The 




Figure 4. Phase diagram of the disordered All class as a func- 
tion of disorder strength W and Fermi level Ef for A = 0.3. 
There is a connected region of quantum spin Hall insulator 
separated from the outer trivial insulator phase by a metallic 
critical region of delocalized states. The width of this region 
in energy is proportional to A (the BIA strength) and is not 
constant as W is varied. 



width of the metallic region in energy appears to decrease 
somewhat over a range of increasing disorder strengths. 
When the BIA term is increased the metallic region per- 
sists over a similar increase in disorder strength, and the 
width of the metallic region in energy also increases by a 
similar amount. 

To summarize, we have shown that the variance of the 
nearest-neighbor level spacing distribution of the entan- 
glement spectrum can serve as a good indicator for re- 
gions of energy which contain delocalized states. We have 
illustrated the nature of the symplectic All class by show- 
ing the existence of wide bands of extended states with 
an energy width which depends on the strength of the 
inversion breaking perturbation. While issues with local 
entanglement due to localized states likely obscure some 
features of the entanglement spectra, a generic feature of 
a flat level-spacing variance is indicative of the Fermi level 
being tuned to a critical region of extended states. This 
observation will be helpful in identifying localized/de- 
localized phase boundaries in disordered fermion systems. 
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